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Abstract 

We address a conjecture introduced by Massoulie (2007), concern- 
ing the large deviations of the stationary measure of bandwidth-sharing 
networks functioning under the Proportional fair allocation. For Marko- 
vian networks, we prove that Proportional fair and an associated re- 
versible allocation are geometrically ergodic and have the same large 
deviations characteristics using Lyapunov functions and martingale ar- 
guments. For monotone networks, we give a more direct proof of the 
same result relying on stochastic comparisons that hold for general ser- 
vice requirement distribution. These results comfort the intuition that 
Proportional fairness is 'close' to allocations of service being insensitive 
to the service time requirement. 

1 Introduction 

Bandwidth-sharing networks describe the evolution of the number of flows 
(or calls) in a communication network where different classes of traffic com- 
pete for the bandwidth. 

They have become a standard modeling tool over the past decades for 
modeling communication networks [U [IT] and have been used in particular 
to represent the flow level dynamics of a wide range of wireline and wireless 
networks [2], generalizing henceforth more traditional voice traffic models [7]. 



In queuing theory, these models boil down to a particular class of pro- 
cessor sharing networks with state-dependent service rates. Assuming that 
class-i flows arrive subject to a Poisson process of intensity Aj and require 
exponentially distributed service times of mean m (the arrival processes of 
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all classes being mutually independent), and in the absence of internal rout- 
ing, the stochastic process X = (Xi, . . . , Xjy) describing the number of flows 
(or calls) in progress in the network is a multi-dimensional birth and death 
process with transition rates: 

q(x,x-ei) = fH<pi(x), 
q(x,x + a) = Xi, 

where x = (x\, . . . , xn) is the number of flows in each class of traffic. The 
service rates of the N traffic classes = encodes the particular- 

ities of the network resulting from the specific topology, technology, radio 
conditions, interference and multi-diversity effects and the packet protocols 
and congestion control mechanisms in use. In a wireline network, the vector 
4>(x) is usually assumed to belong for all x to a polyhedron describing the 
capacities constraints of each link that are used by the different routes while 
for wireless networks, (j)(x) generally belongs to a more complicated closed 
convex set containing corresponding to the achievable rates. In many situ- 
ations like TDMA or CDMA data networks, the capacity set can be assumed 
to be convex but this is not necessarily the case for decentralized schemes 
like for the 802.11 WLAN. However, it has been shown that the capacity 
sets in that case are log-convex (see [9]). 

Some specific bandwidth allocations have received a lot of attention in 
recent years. These include the max-min fair allocation and the proportional 
fair allocation (PF) that maximize a log-utility function. More generally, Mo 
and Walrand introduce the following family of utility functions: 

N i_ Q 

[/ ft (j;,i|) = Vi^ , a£(l,oo), 

^— ' 1 — a 

i=l 
N 

U 1 ^,!]) = ^2xilog(rn), a = l, 
i=i 

including as special cases the proportional fair allocation and the max-min 
fair allocation for a tending oo. On the other hand, the balanced fairness 
allocation (BF) was defined in [3] as the allocation ensuring the reversibility 
of the Markov process X and maximizing the probability of the network to 
be empty (among the 'reversible' allocations). 
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More generally, the selection of a specific bandwidth allocation is moti- 
vated by several properties of the resulting process X. Among those prop- 
erties, the following ones are of particular interest: 

1. maximal stability: one expects the bandwidth sharing mechanism to 
stabilize the system whenever it can be stabilized, 

2. decentralized protocol: the bandwidth allocation can be implemented 
in the network using decentralized schemes, 

3. robustness: when changing the traffic condition^ an d in particular 
the service time distribution (but keeping the mean flow size fixed), 
one could expect the stationary measure of X to remain the same, in 
which case the system and the bandwidth allocation are said to be 
insensitive, 

4. relaxed robustness : the large deviations characteristics do not depend 
on the service time distribution, except for its mean and coincide with 
the large deviations of the most efficient insensitive allocation. 

The PF allocation satisfies properties (1) and (2) but fails to satisfy 
property (3) on general topologies. On the other hand, the BF allocation 
satisfies (1) and (3) but it is not known whether it satisfies (2) [2]. On the 
one hand, it was shown in |10| that an appropriate modification of the pro- 
portional fair allocation, called modified proportional fair allocation, (mPF) 
coinciding asymptotically (point-wise) with PF, has the same large devia- 
tions characteristics as BF. On the other hand, it has been recently proven 
that an insensitive allocation being maximal stable, is asymptotically equiv- 
alent (point-wise) to PF j!5j . However, it remained an open problem to 
prove that the large deviation characteristics of the stationary measure of 
the PF allocation itself coincide with those of mPF, and BF as it was con- 
jectured in [TO] . 




Our contribution is the following. We first recall the principles of the 
Freidlin and Wentzell theory for birth and death processes on 7j N with rates 
being Lipschitz and with bounded logarithms. As underlined for instance 
in |14j . it is very demanding to extend these results to processes on state 
spaces with boundaries since the technical conditions of the classical theory 

1 When the size distributions of the flows are not exponentially distributed, the process 
X is not Markov by itself anymore and the dynamics have to be defined using the residual 
service time. 
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are never fullfilled. (Remark that boundness of the logarithm is never veri- 
fied for birth and death processes living in the orthant, while the Lipschitz 
assumption is not verified in our case for several network topologies). 

A main contribution of the present article consists in overcoming these 
difficulties for the specific processes we are studying. We first show that for 
Poisson arrivals and exponentially distributed flow sizes, the stationary dis- 
tribution tt pf of the number of flows associated with PF and the stationary 
distribution -jr mPF associated with the mPF allocation have the same large 
deviations characteristics. More precisely: 

Theorem 1. For all x G Tl^l : 



In the particular case that x G the bound can be improved to 0(n 1 ). 

This is achieved by first proving the geometric ergodicity of both the 
mPF allocation and the PF allocation. For that purpose, we exhibit ap- 
propriate Lyapunov functions, relying on some structural results of PF de- 
scribed in [10] . This then allows us to use simple martingale arguments. 

Finally, for monotone networks and generally distributed flow size, we 
give a more direct proof establishing that the large deviations characteristics 
are actually insensitive to the service time distribution. This shows that the 
proportional fair allocations indeed satisfy properties (1) and (2) and (4) at 
least on monotone topologies. 

The rest of the paper is organized as follows. In Section [5J we intro- 
duce the allocations functions used in the sequel as well as several of their 
properties. In Sectional we show that the rate function of the proportional 
fair allocation coincides for Markovian dynamics to the rate function of Bal- 
anced fairness. In Section [H we prove property (4) for monotone networks 
using stochastic comparisons. 

2 Properties of bandwidth sharing allocations 

Notations We define here a few notations that we need in the sequel. 
M+ denotes the set of non-negative real numbers. The bandwidth allocation 
vector is denoted <f>{x) = (^i(x))i=i...7v '■ For any vector v in M w , and function 




Ve > 0. 
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/ : R — > R, we denote f(v) the vector (/(«j))i=i..jv- Similarly, we also use 
the notation a" for ]X =1 a T • The usual scalar product in R^ between u and 
v is denoted {u,v}. We use 1 1 • | | p to denote the Z p -norm, but we reserve | • | to 
denote the Zi-norm: \v\ = Yli=i \ v i\- For x i U we a ^ so use the notation 

x < y to denote the partial order Xi < yi for all i = 1 . . . N. Finally, given a 
set C and function / : R — )■ R, we shall denote /(C) the set {/(??) : i) 6 C}. 

2.1 Gradient allocations 

In the characterization of the large deviations of the stationary regime, we 
shall rely on two properties of the allocation function (j) playing a crucial 
role and being intrinsically related: being a gradient allocation or a discrete 
gradient allocation. 

Definition 1. A bandwidth allocation is called gradient if there exists a 
function P : M. N — > R (that we call a continuous potential) such that 

log(0(x)) = -VP(i), VxeR^\{0}. 

Recall that the proportional fair allocation is defined by the optimization 
problem associated to a capacity set C as follows: 

(f) PF (x) = argmaxC/ 1 (r/,x) = argmax(a;, log^)}. 

Following Massoulie [ID], observe that the proportional fair bandwidth 
allocation is gradient. Indeed, let 8^ the support function of a bounded 
convex set A i.e.: 

$a( x ) = max i x i ij)- 

Proposition 1. Assume that the set C is log-convex, i.e., the set log(C) is 
convex, then 

log(^ PF (x)) = V8* log(c) (x), V* € \ {0}. 

Proof. 

The function 8* is sub-differentiable because it is convex and finite (see 
|13| ) for all x £ R^. The unicity of the sub- gradient comes from the strict 
concavity of the log function and implies the differentiability. 

□ 

In the sequel, we always assume that C is convex and contains the set 
{v '■ Vi — c } f° r some c > 0. Of course, this is not a restriction for 
applications. We denote by P PF = 8f , C \. 
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2.2 Reversible allocations 

Definition 2. A bandwidth allocation is called reversible or discrete gradi- 
ent, if there exists a function P : H. — > M. (that we call a discrete potential) 
such that for all x € N N , i = 1, . . . , N: 

log(^(x)) = -DP(s) = P(x) - P(x - a). 

In the latter case, the stationary measure of the process is easily de- 
scribed: 

Proposition 2. The process X associated with a reversible allocation with 
discrete potential P is reversible ( in the usual sense ) and its stationary mea- 
sure (when it exists) is: 

ir(x) = C\ x exp(P(x)). 
Assume further that -P(nx) —> j(x) as n — > oo. Then: 

1 N 
lim -logins) = -(7(2;) - V" Xi log(Aj)) . 

n— >oa n * — ' 

j=l 

A particular role in what follows shall be played by two reversible alloca- 
tions 'close' to PF: the modified PF allocation (mPF) and the BF allocation 
which we now define mathematically. 

Definition 3 (Modified proportional fair and Balanced fair allocations). 

log(<f> mPF (x)) = -VP PF (x), 
with P PF the continuous potential associated with PF, and 

log(^ BF (x)) = -DP(x), 
where the potential P is recursively defined by: 

P(x) = 0, 

P{x) = max{a > : P(x - a) - a € log(C)}. 
We have the following large deviations results for the BF allocation: 
Proposition 3. (Massoulie, 2007) 

lim -log7r BF (nx) = P mPF {x). 

n— >oa fl 
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The main challenge of the subsequent analysis is to provide similar re- 
sults for the Proportional fair allocation itself and hence verify that both 
allocations share the same large deviations characteristics. The main diffi- 
culty of this program is that the stationary measure of a network under the 
proportional fair allocation does not have a closed form in general (except 
in the very particular cases of symmetric hypergrids). 

2.3 Stability 

We state two stability results, one for the original stochastic system, and one 
for the deterministic analogue Both results are well known for Markovian 
bandwidth sharing networks in pQ. Massoulie extended the stability results 
to general service time distributions in |1U| . 

Theorem 2 ([10]). The stability set of the process X associated with the 
PF allocation is contained in, and contains the interior of the set S = C. 

Remark 2.1. For Markov processes, a direct proof follows from taking x i-4 

YliLi X" as a Lyapunov function. 

Define now the deterministic dynamical system (DS) XQ by 



It was proved [8] (see also [10]) that the fluid limits of the process X are 
solutions of this dynamical system at all regular points. 

2.4 Examples 

We illustrate the above description on a few toy examples. 

Example 1 (Single link). Consider first a single link shared by N classes. 
This corresponds to choosing C = < 1}, where the {//?} correspond to 

the mean flow sizes. In this case, the proportional fair allocation coincides 
with the balanced fairness allocation. Its service rates are given by: 

4>i(x) = Mirr- 
\x\ 

Note however that the continuous and the discrete potential functions asso- 
ciated with the allocation (ft do not coincide: 



(Xi - <f>i(x))l Xi>0 




x(0) =x . 
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while the discrete potential is given by: 



P(x) = log ( ( 1X1 

\\Xl, ...,X N 

One can however prove using the Stirling formula that 

5 *og(C) 



— (P(nx) — 5j 



An extension of these formulae can be obtained for hypergrids topologies 

Example 2 (Tree network). A 2-level tree network with N routes is defined 
by the following polyhedron for the capacity set: 

P = { n sR^:r H <c i , < 1 }, 

i 

for some constants {q}. An expression of the stationary measure for the 
balanced fairness allocation can be found in In general (i.e., when the 
tree is not degenerate), the proportional fair allocation does not coincide with 
the balanced fair allocation and does not raise a reversible process. Hence 
its stationary distribution can be obtained only numerically. For the case 
N = 2, and max(ci,C2) < 1, the proportional fair allocation is: 

<j>\{x) = 1 - (p 2 (x), 

X2 \ 



(p2 (x) = min(c 2 , 



X\ + X2 

An important property of the proportional fair allocation in that case is 
its monotonicity (see Section 4 for a definition) |5j/; 



^f F {x) . 4>f F {y) , w . 

> , Xj < yj, Vj. 



Xi y% 

Example 3 (Wireless network). In J2j, the rate region C of a 2- station 
network functioning under the 802.1 le protocol is studied. The rate region 
is the set of achievable throughput vectors at the fine time scale (packet level). 
Their findings show that the rate region, ( which exact expression depends in 
a complicated manner of the probabilities of transmitting) is generally not 
convex but is however log-convex. 
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3 Large deviations for Markovian dynamics 



3.1 The Freidlin-Wentzell theory for birth and death pro- 
cesses with smooth transitions 

In this Section, we remind the concepts introduced by Freidlin and Wentzell 
[6] allowing to get a grasp on the sample-paths large deviations of stochastic 
differential equations with small noise. These results were often applied to 
diffusions processes but can equally be applied (as indicated in [6]) to some 
jump processes. We show here how these classical results can be used to 
prove a large deviation principle for birth and death processes with "gradi- 
ent" rates on "L N . 

Let Y n be a multi-dimensional birth and death process with transition 
rates: 

/x x e;\ fx\ 
<?(-, J = n<pi[-), 

(X x e;\ 

q -,-H = nXi. 

\n n n / 

We suppose additionally that the death rates are O-homogeneous (i.e. <fi(az) = 
<p{z), Va > 0). 

Define the logaritmic moment-generating of the increment of the process 
for z € M. N /n,y G R N by 

H n (z,y) = j t E*[eM(y,Y n (t)-z))]. (3) 
Using the structure of the generator and the 0- homogeneity of the rates: 

H n (z, y)=n(^2 ^ mln " 1) + &(*)(e- w/n - 1)1 • (4) 

We thus obtain that H n (z,y/n) = nH(z,y/n) with 

H(z,y) = ( e y + (e~y - l,<f>(z)). (5) 

Define now the Fenchel-Legendre L transform of H: 

L(x,y)= sup(y,e)-H(x,e), (6) 

eeR N 
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and define the action functional S : C([0, T]) — >• K as follows: 



Sr(r) = / L(r s ,r s )ds. 



(7) 



./o 



Now define the quasi-potential V by: 



V{x) 



mm 

r,T,r(0)=0,r(T)=x 



S T (r,r). 



(8) 



Assume that the function L is such that: 



sup 

\x— x\ 



L(x,y) - L(x,y)\ 
1 + L(x,y) 



0. 



The original results of Freidlin and Wentzell are stated for Lipschitz- 
continuous transitions (extended to W N ): 

Theorem 3 (Freidlin and Wentzell). Assume that for each i, the junction 
log^)j(-) is bounded and Lipschitz continuous. Let ir n the stationary measure 
of the birth and death processes X n . Then, a large deviations principle holds 
for the family of probabilities ir n (-) with rate function V(-). 

3.1.1 Further characterization of the potential 

In the case of general multi-dimensional birth and death processes, it is dif- 
ficult to solve the variational problem from which the potential V is defined. 
However, assuming that the rates are gradient, as defined in the previous 
Section, greatly simplifies the expression of the potential. 

Proposition 4. Assume that the rates (allocation) <f> is gradient with con- 
tinuous potential P. Then the quasi-potential of the birth and death process 
is equal to R(x) = — (log(A),x) + P(x). 



The first simplification comes from the time homogeneity of the process. 
Using the Euler-Lagrange principle, combined with the time homogeneity 
(which implies that V X L = 0) the minimizing path satisfies the Beltrami 
identity: 



Proof. 



We first show that: 




L-(y,V y L} = 0. 
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Now, observe that if 9 y = argmax(y, 9) — H(x, 9), we obtain that V y L = 9 y 
and L — (yV y L) = H(x,9 y ) = 0. Solving the last equation leads to: 

e± = 2 log (^~)' 

which allows us to conclude since 

L(x,x) = (x,V y L(x,x)) = (x,0&). 
The expression for the quasi-potential then follows from integration. 

□ 

We insist on the facts that this theory allows to prove the large devi- 
ation rate of the stationary measure of birth and death processes only for 
smooth allocations and that it becomes a very technical issue to weaken this 
assumption. Also it does not provide a rate of convergence in general. 

We address these issues in the next Section: we prove directly the large 
deviation principle using martingale arguments by taking advantage of the 
geometric ergodicity of the process. 

3.2 Large-deviations results for general gradient allocations 

In this Section, we call X the process associated with the PF allocation, 
denoted 4> an d corresponding to the continuous gradient of the potential 
P PF . We further denote by X the process associated to the allocation <fi 
itself corresponding to the discrete gradient of P PF , i.e.: 

logOK*)) = VP PF (x), (9) 
log(0(x)) = -nP PF (x) = (P PF (x)-P PF (x-ei)) l=1 ... N . (10) 

Structural properties of the PF allocation 

From the structural representation of PF, we know that 4> is a perturbation 
of the original rates <j) in the sense that (see Lemma 9 in [10] ): 

|log(&(s))-log(&(x))| = \log( ( / )l (x))-(P PF (x)-P PF (x-e i ))\ < -, Vt,x. 

Xi 

(11) 
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Geometric ergodicity 

We first consider the dynamics corresponding to the Proportional fairness 
allocation. 

Proposition 5 (Geometric ergodicity of X). Suppose A G int(C), then there 
exists a constant K such that for \x\ > K : 

AG(x) < ~lG(x). 

Hence: 

P?{x) - 7T PF ( X ) < K ie - K2t , 

for some constants K\,Ki > 0, and ir PF the stationary distribution of X. 
Proof. 

The proof has two steps. We first construct a Lyapunov function with 
bounded drift. We then use this Lyapunov function to construct a new one 
verifying a geometric drift inequality. 



First step: let F{x) = (J^Li y - ) 1 ^ 2 - Observe that F is a norm in M. N 
(hence is positive, 1— homogeneous, and diverges to infinity when |x| — > oo). 
Furthermore, it is C 2 for all x ^ and: 

VF(x) = -^rr . 

Hence, there exists K > such that for \x\ > K 

1 ,d 2 F(z) l 
|— I sup 1 < e. 

\X\ z=x,x±ei O Xi 

Using that F is 1-homogeneous, this leads for |x| > K to: 

N 

AF(x) = ^A i (F(x + e i )-F(x)) + ^(x)(F(x-e l )-F(x)), 

i=l 

< (X-(/)(x),VF{x)} +e, 
= + 
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Using the definition of Proportional fairness (as the allocation maximizing 
U (-,x) and the strict concavity of the log function), for all r\ € C, 

~^Vi Vi 
which implies that that there exists 7 such that: 

(A — <p(x) } i— 7T-) < -7, 

which in turn implies (together with the fact that F is norm-like) that there 
exists 7 such that: 

AF(x) < -7, Vx, |x| > K. 

Remark also that 

\AF(x)\ < C, Vx. 

Second step: we now calculate the drift of G(x) = exp(5F(x)): 

N 



AG(x) = ^A i (e 5F ( x+ei )-e F ^)) + ^(x)(e 5F ( a; - e ')-e 5F ^)), 
i=l 

N 

= G(x) ^ X^e^x+^-FW) - 1) + <t>i(x)(e s( - F( - x - e ^- F ^ - 1)), 

N 

< G(x)(5AF(x) + Ci<5 2 ^exp(c|AF(x)|) + 0(<5 3 )), 

< G(x)(-C 2 <5 + C 3 <5 2 ), 

< -7G(x). 

This implies that (see j!2j): 

|P°(X(t) = x) -vr PF (x)| < ^ 3 exp-^ 4t . 

□ 

We now deal with the dynamics of the modified proportional fair alloca- 
tion. 
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Proposition 6 (Geometric ergodicity of X). Suppose A G int(C), then there 
exists a constant K such that for \x\ > K : 

AG(x) < -lG(x). 

Hence, 



P?(X)-TT PF (X) 



-K 4 t 



for some constants K%,K^ > 0, ir mPF being the stationary distribution of 
X. 

Proof. 

We proceed as previously. The only difference consists in proving that 
the perturbations of the rates are small enough to be negligible in the drift 
calculations. Let F(x) = {^=1 jt) 1/2 - For M > K, using the bounds on 
the difference of 4> and (ft, we obtain that: 

\(fti(x) - 4>i(x)\ < cq\ log(0i(x)) - \og(4>i(x))\ < — . 

X{ 

Hence: 

N 

AF(x) = ^A i (F(x + e i )-F(x)) + ^(x)(F(x-e i )-i ? (x)), 

i=l 

< (\-j>(x),VF(x)) + e, 

< (X-<ft(x),-^-) + ( ( ft(x)-4>(x),-^-) + e, 

F(x)X F(x)\ 

< -7+(0( s )-0(x) _ ^— ), 

F(x)\ 

c 

< -7 + rT- 

\x\ 

which gives (together with the fact that F is norm-like) that there exists 7 
such that: 

AF(x) < -7, Va?, |x| > K. 

Remark also that 

\AF(x)\ < C. 

We can conclude following exactly as in the previous proof. □ 

Before proving the main result of this section, we establish a useful 
lemma using the geometric ergodicity of the processes. 
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Lemma 3.1. There exists C > such that for t n = Cn, we have 

(X tn > nx) 



log 



TT PF ({nx}^) 



< C\ exp{— C2 n} 



for C\,C 2 > 0, where ir PF \{k}^) = X]j>fe 'K PF (j). The analogue inequality 
is true for process X. 

Proof. 

We need to bound from below the stationary probabilities of the process. 
This is easy since the rates are bounded. Let us define a process with the 
same arrival rates and the service rates equal to 4>, the maximum of 4>i for 
all coordinates. It is clear that 



for some constant K\ > 0. 

We have 

w0 (X tn > nx) 



log- 



vr PF ({nx}t) 



< 



t n > nX ) 



7r PF ({nx}t) 
^{X tn > nx) - Ti PF ({nxY) 



Tr PF ({nx}T) 



< 



< 



K 2 e~ K3tn 
7r PF ({nx}t) 



K 2 e 



-K 3 tr, 



Mi) 



K 2 (<t>\ 
— exp{-K 3 t n + log 



where first inequality comes by log(x) < |1 — x\ for all x > 0, and the second 
one by the ergodicity of X. We take C > such that K3 C > log (^). 
Note that the same argument functions for process X taking its respective 
transition rates and stationary distribution. □ 



Change of measure and control of the martingale 

To relate the distribution of X and X, we recall the Proposition B.6 513 of 
Schwartz and Weiss [14 



Proposition 7. Let Y£ fori = 1,2 two multidimensional birth-death process 
in with step directions ej, bounded step rates {q){x)} x ^N and law P i . 
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Assume that for all x and j, q^{x) = if and only if q 2 (x) = 0. Then we 
can relate the distributions of the processes by dP 2 = MtdP 1 where 




Here {N J } denote a family of counting processes describing the jumps of the 
processes in the j-th direction and Mf is a cadlag martingale. 

In our case, (fii(x), 4>i(x) are both positive whenever Xi > 0, and oth- 
erwise. All the rates of X and X are bounded so we meet the conditions of 
the previous Proposition and we can write that dP(u) = AltdP(u), with 

M t = ex P (/EW>o } log(||^)^ 
-J hx s >o}(M x s)-M x s))ds), 

N 

= U M i> 

and each of the M 3 is itself a martingale. 

Remark 3.1. An important observation for the following is that the count- 
ing processes N 3 can be seen as thinning ( according to X s ) of some Poisson 
processes N J which are all independent. 

We denote by M s ^ the last expression with the integral running from s 
to t, so M t = M , t . 

The change of measure formula is then: 

E[l { x t=nx} ]=E[M t l {Xt=nx} ]. (12) 
We can now prove our main result. 

Proof of Theorem 1: 

We consider first the case where Xi > for all i = 1 . . . , N. 
Define the sequence of stopping times: 

n = o, 

Tj = inf{i > Ti^\,Xf > nx}, for i even, i>2, 
Ti = inf{t > Tj_!,X t < nx}, for i odd, i > 3 
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Observe that if Xt > nx, then necessarily, there exists k even (a.s. finite) 
such that Tfc < t < Tfc + i. Using the Markov property and the martingale 
property, we then have 



E °( M 0A{X t >nx}) ~ E ° Yl 1 ^<*< T -fc+i M * 1 {X t >nx} 

k even 

= ^2 EOl T k <t<T k+1 M 0tTk M Tktt l { x t > nx} , 

k even 

= Yl S °(^( 1 ^<*<r fc+1 ^0,r fc M Tfc , t l { x t > ra} l^r fe )) 
k even 

= E E °( M Q^ k E (U k <t<T k+1 M Tktt l {Xt > nx} \X Tk )), 

k even 

with {J^k} being the natural filtration of process X. 
We define 

9k = supE{l Tk < t<Tk+1 M Tk:t l { x t>nx} \X Tk = y), 
v 

so 

E °(M t l {Xt > nx} )< £ 9kE°(M , Tk )= £ g k . 

k even k even 

Realize that on {X t > nx} D {r^. < t < t^+i} (k even) we have {X s > 
nx : s E [Tfc,i]}. Recall that x is such that Xi > for all i. Using the 
assumption on </> and 0, we get that asymptotically in n: 

E ° (K<t<r k+1 \ Xt > nx} ex P ( £ £ i {Xs _ >0} lo s(|j|^y W 

-/ T, 1 {x.->o}fa( X ')-M x *))d* 

Jr k j=l 

N c 

< E °{K<t<r k+1 e X p{^ + t)}\x t >n X })i 

3=1 



< E ° (lr k <t<r k+1 eX P {^(7V t + t)}l {Xt ^ 



CiiV 



n 



L {Xt>nx} 



where Nt is a Poisson process with parameter A equal as the maximum of 
the parameters of the Poisson processes {Nf}j =1 . 
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Summing these inequalities and using Holder's inequality 



n Cl Nt n ~ P- 1 - CipN 1 

E°(M t l { x t > nx} ) < exp{^— }¥°{X t > exp{At[e^ - 1]}?. 



We now choose i n = C^n such the result of Lemma 13.11 holds, and the 
sequence p n = n to obtain 

1 ^ pCVv ^ , . C X NC 2 , AC^exp^AQ - 1) 

-logP u (X tn > nx) < 1 

n n n 

+ ^logP (l tn >nx), 

-log(-K PF ({nx}T)) + ^exp(-C 4 n) < 0{n~ l ) + ( - log(7r mPF ({nxP 
n n \ n J n 

, C 3 (n-1) ^ 
H 5 exp( — C4 nj . 

So we have 

-logfTr^anx^)) < ilog(vr mPF ({?ix} t )) + 0(n- 1 ). 
n n 

We also have the converse of last inequality using the same arguments in- 
terchanging X and X. 

Assume now that Xj > for i £ W, while x% = for % £ S, where 
U, S C {1, iV}. If classes of W are independent of classes <S, the result is 
obvious. Hence suppose that there exists a coupling between classes of U 
and S. We shall make use of the following fact: 

Lemma 3.2. Given the definition of Proportional fairness, if class i € S 
and class j GU are coupled, then there exist K > such that if yi < ne and 
Uj > n for all j GlA, then: 

Mv)-Hv)\ <Ke, 

Proof. 

If at least two classes are coupled, there exists a Lagrange multiplier 
a > and some positive constants c% > and c,- (with at least one j such 
that Cj > 0) such that 

dU l (y,r)i) Xi 

i = — = aci, 

or\i rji 
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drij Vj 



acj. 



Combined with the fact that there exists c > such that < c < Ej=i Vj — 
C, we obtain that 

<Mv)<c- 



Vi + EjeW> Vj ~ Vi + Ej-eMi Vi ' 

where Ui CW. Hence if yi < ne and j/j > n for all j Gti 

~ f n ~ 

My) < c-— < Ke. 

1 + n 

Using the control on the modified proportional fair allocation for Xi > 1: 

My) < ex P (—)My) < ^My) < e-ke. 

□ 

Now we do a finer classification of indexes in S. Let {e n } be a sequence 
such that e n goes to as n grows, and {t n = C2T1}. For each n, we define 
the sets 

S n = {i G S : Xi(t n ) > ne n }, S n = {ieS: Xi(t n ) < ne n }. 
We look at the set of events: 

A = {Xi(t n ) > nxi,i G Ur, Xi(t n ) > ne n ,i G <S n , Xi(t n ) < ne n ,i G S n }. 
Recall that 

|los W £1 ' 

and remark that on A G A, using the previous lemma, the process counting 
the number of downwards jumps in direction j G S which are not common 
for X and X is dominated by a Poisson process N 6n of intensity C§e n inde- 
pendent of N 3 , j G S . 
Hence: 

E°l Tk < t<Tk+1 l A M Tk j 

< E°(e W (J2—(Ni + t)+ — (Ni+t)+ J2PsNt- + (C 6 e n )t])l Th < t<Tk+1 l A ). 
W n je s n n€n jeS n 
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and, 

E°(M t l A ) < £° ( exp (V — (Ni + t) + £ -^(iV/ + t) + ]T + (C 6 e n )t]) U 

< E o ( exp (#UCi(N t + t) + ^^(N t + t) + #S n [C 5 N^ + C 6 e n t]) U 

V \ Til Tl/dfi ' 

<^(exp { ^(iV 4 + t)})-^(e X p { ^(iV t + i )})^ 

V n /V ne n / 

E°(exp{ P N{C 5 Nt + C 6 e„t]})^P°(Xt > nxfr , (13) 
where we used Holder's inequality for 1 = |^ + ^ + + 2^1 . 

We now bound the speed of convergence, in the large deviations scale, 
of each of the three first factors on the right-hand side of inequality (|13p. 
Choosing t n = C 2 n, p n = Cjlogn, e n = n~2 i we have: 

-log£ u exp{ (iV in +t n )| = -log exp{ }exp{Ai n e » - 1 

n V n J n \ n 



= #ZYCiC 2 + m rte _ ^ 
n Nnp n 

_ mAC\C 2 #U\C 2 C1C7N log n 

~ n + JVCVlogra^ ^ 
= 0(n" 1 ). 

We do similarly for the second factor: 
-log£ u exp{ W„+*n)} = -log(exp{ }exp{At n [e^« -lj}^ 



n e n n / n \ e n n 



1 e e ™ n — 1 

e n n Nnp n 

e e " n —1 

e n n NC 7 log n L J 

©((ne,,)- 1 ) =0(n"2). 
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Finally, for the third factor we have: 



1 / \ 1 / #5 

log E°(exp{N Pn [C 5 N% + C 6 e n t n ]}J Npn = -log ( exp{C 6 #S n e n t„} exp{e n t n [e^^ _ 



n \ / n 

n Np n n 
= C,C 2 #S n e n + clft^ 0507 " 10 ^ ~ $ 
rr ,o f i C 2 #S n e n [n c ^ N - 1] 

- CtCrfs^ + C7iVlogn 

= 0( en ) + 0(^l) = 0(n-i+ e ). 
logra 

where we chose Cj > small enough such that CjC^N < e, for any given 
e > 0. 

Then along the same lines as in the case x has all entries positive, it 
follows 

- log(n PF ({nx}T)) < - log(vr mPF ({nx} t )) + 0{n~^ +e ) Ve > 0. 
n n 

□ 



4 Large deviations for monotone networks with 
general service time distributions 

In this Section, we consider a processor sharing network (i.e. a set of pro- 
cessor sharing nodes) with a proportionally fair bandwidth allocation. This 

means that a flow is served at node i with speed ^r^j^ • We need the 
notion of strong monotonicity introduced for instance in [5]. 

Definition 4. (j) is strongly decreasing if the function tpi defined by ipi{x) = 
i^Sl if Xi > and otherwise, is decreasing in Xj for all j ^ i. 

We assume in this Section that the network is monotone, i.e. that <p is 
strongly increasing. This is verified on all tree topologies for instance and 
for many wireless networks instances. In the proof of the following theorem, 
we use the monotony of the network to obtain stochastic comparisons. We 
recall the following Proposition: 
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Proposition 8. Let Xt and Xt two processes associated with the allocations 
4> and d>. Suppose that <j) is strongly decreasing, and that ^iM > Mg) t or 
each x. Then for all t and for all service time distribution: 

Xt -<st x t . 

Proof. 

Remark that the processes Xt and Xt are not Markov in general. We can 
however easily construct a sample-path comparison between the processes. 
Starting with ordered initial configuration of flows x -< y, the rates at which 
each flow is served in the network verify: 

4>i(x) > My) > My) 

%i Vi Vi 

Hence, a coupling can be constructed such that X t -< X t almost surely. 

□ 

Recall that we define P PF as P PF (x) = max^ g i og ( C )(j],i;}, and that 
log(&(x)) = VP PF (x). 

Now define the discrete potential functions ^ and by the recursive 
formula: 

^(x) = max (\S>(x — e^) — log(0j(x))), 

i=l...N 

^i(x) = min {^_(x — ei) — log(4>i(x))} . 

i=l...N 

We call supPF and infPF the reversible allocations associated with the dis- 
crete potentials and and define the rate function R(x) = (log(A),x) — 
P{x). We can now state the main result of this Section: 

Theorem 4. Suppose the network monotone, i.e. 4> is strongly monotone. 
If X is in the interior of the capacity setC, then the allocations Proportional 
fairness, supPF and infPF as well as Balanced fairness are stable and admit 
the same large deviation characteristics with rate function R. 
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Proof. 

Given the definitions of the discrete potentials and observe that 
there exists two paths V{x) and V_(x) from to x such that: 

= - Yl iog(^ fc (^)), 

k:(i k ,z k )eV(x) 
k-{i k ,z k )eV{x) 

where with a slight abuse of notations the indexes i{z) correspond to the 
indexes defined by the specific paths V{x) and V_(x). 

Let X and X_ the process associated with the balance function ^ and 
\£. Recall now that P is the potential associated with the proportional fair 
allocation, and since we can define this potential up to an additive constant, 
let us choose it such that P(0) = 0. So, we can write that for any path 
V(x) going from 0_to x, P(x) = P(x) - P(0) = Ek:(i k ,z k )eV( x ) V ik P{z k ). 
Choosing V{x) = V(x): 

\*(nx) - P(nx)\ = | Vi k P(z k )-log(& k (z k )\, 

k:(i k ,z k )£V(nx) 

< Yl \D lk P{z k )-\og{^ k {z k )V (14) 

k:(i k ,z k )£V(nx) 

Recalling that for Xi > 1: 

\P(x) - P(x - et) - log(^(x))| < — , 

Xi 



we can bound (|M[) by: 



\V(nx) - P(nx)\ < ^ ~~i 

(i k ,z k )eV(nx) %h 



n 

< Cj2^<C'\og(n). (15) 



i=i 



Define the invariant measures (not necessarily stationary) 7r and n of the 
processes X -< X -< X_. Since we have constructed these processes from a 
discrete potential, we have that (see Proposition [2]): 

lf(x) = \ x exp(-$(x)), 
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7r(x) = X x exp(-*(x)). 

Now observe that: 

1 _ l_ 

— log(7r(nx)) = (log(A),x) ^>(nx), 

n n 

— log(7r(nx)) = (log(A),x) — — ^_(nx). 
Using (|15|) we obtain that 

— log(7f(rex)) = — (log(A), nx) P(nx) + 0(n _1 ), 

Th Th Th 

— logMnx)) = —(log(X),nx) — — P(nx) + 0(n~ l ), 
n n n 

and hence, by the 1-homogeneity of R, 

lim — log(7f(nx)) = lim — log(7r(nx)) = R(x). 

n— >oo n n— >oo n 

Assume A is in the interior of the capacity set C. Using the definition of 
P, this implies that there exists a > such that 

R(x) = (log(A),x) -P(x) = (log(A)-log(0(z)),x) < -a\x\, 

which implies that the invariant measure W and tt_ are summable. This 
hence proves that the stationary distributions of both processes X and X_ 
are well defined for Markovian dynamics while the reversibility condition 
(i.e. the balance property of the service rates) implies the insensitivity of 
the stationary distribution to the service time distribution J3l [16] . Hence, 
for any service time distribution, the networks with allocation infPF and 
supPF admit the stationary distribution C\T and C^TL where C\ and C2 are 
normalizing constants. Now using the assumption of monotonicity of the 
network, we obtain that: 

P(\X\ >n)> P(\X\ >n)> P(\X\ > n), (16) 

and by Cramer's theorem we conclude: 

lim -logP(\X\ > n) = lim - log P(\X\ = n). 

n— >oo n n—too n 

It has further been proven in [10] that mPF and BF have R(x) as rate 
functions. 

□ 



24 



5 Conclusion 



We proved that the stationary measure of the number of flows in progress 
in a bandwidth sharing network functioning under the proportional fair al- 
location shares the same large deviations characteristics with the stationary 
meausure of the number of flows in progress of the same network under the 
balanced fair allocation. This formalizes the idea that in long excursions 
proportional fair allocation behaves similarly to the most efficient insensi- 
tive allocation. 
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